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Abstract. Quarkonia (3/tp, ip', xc, T) production provides a sensitive probe 
^ ' of gluon distributions and their modification in nuclei; and is a leading probe 

1} . of the hot-dense (deconfined) matter created in high-energy collisions of heavy 

ions. We will discuss the physics of quarkonia production in the context of 
recent p + p measurements at PHENIX. We next discuss Cold-Nuclear Matter 
(CNM) effects as seen in our measurements in d + Au collisions - both for 
intrinsic physics such as gluon saturation and final-state dissociation, and as 
a baseline for studies in nucleus- nucleus collisions. Then we review the latest 
nucleus-nucleus results in the light of the expected CNM effects, and discuss 
■^j- ' two leading scenarios for the observed suppression patterns. Finally we show 

the latest data from PHENIX, including new d + Au data from the 2007-2008 
run; and then look into the future. 

^■O Keywords: quarkonia, heavy-quarks, quark-gluon plasma, gluon saturation, 

cold nuclear matter. 
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1. Introduction 

We discuss our present understanding of Quarkonia (J/ip, ip' , XCi T) based on 
the measurements by PHENIX at RHIC. We discuss 1) production, 2) cold nuclear 
matter (CNM) effects, 3) the effect of the Quark Gluon Plasma (QGP), and then 
comment on prospects for the future as RHIC luminosities increase and detector 
upgrades are installed. As shown in Figure[TJ the numbers of J/ip obtained in recent 
runs has increased dramatically, with over 70,000 in the just completed d + Au run. 



2. How are Quarkonia Produced 



Quarkonia are produced primarily via gluon-fusion, but it has proven difficult for 
theoretical predictions to reproduce both the cross section and the polarization of 
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Fig. 1. Approximate Number of J ftps per year for different types of collisions at 
PHENIX. Close symbols are for dimuons at forward rapidity, and open symbols are 
for dielectrons at mid rapidity. 



the J/tp- The configuration of the initially produced cc state remains unclear, and 
casts uncertainty on what CNM effects it will experience in nuclei. NRQCD models 
produce a cc in a color-octet state and are able to reproduce the cross section, but 
predict large transverse polarization at large pr - unlike the data from E866/NuSea[ 
[T] and CDF [[^ which show only small longitudinal polarization. However, a recent 
color-singlet model[|3] claims good agreement for both cross section and polariza- 
tion. 

Another complication in quarkonia production, particularly for the J/tp, is that 
about ~ 40% of the J ftps come from decays of higher mass resonances, namely the 
ip' and xc- Until recently, these fractions have been inferred from measurements at 
other energies. Now PHENIX has started to quantify these itself with initial results 
indicating 8.6 ± 2.5% from the ip 1 and < 42% from the \c- Another PHENIX 
measurement [ [5] shows that 4^% of the J ftps come from decays of B-mesons, a 
contribution which is strongest at larger px ■ 



3. What Cold Nuclear Matter (CNM) Effects are Important 

For Quarkonia produced in nuclei, e.g. in p+A or d+A collisions, several interesting 
effects - usually called cold nuclear matter (CNM) effects, can occur. These include 
modifications of the initial gluon density either according to traditional nuclear 
shadowing models [ [5] that involve fits to deep- inelastic scattering and other 
data, or gluon saturation models [0]. In addition the initial-state projectile gluon 
may lose energy before it interacts to form a J/tp. Both of these effects can cause 
suppression of the produced J/ips per nucleon-nucleon collision at large rapidity 
(or small x) relative to that observed in p+p collisions. Finally, the J ftps can be 
suppressed by dissociation of the cc by the nuclear medium in the final state. 

A new analysis of the 2003 PHENIX d + Au data, along with the new 2005 
baseline p+p data have been put together to produce new nuclear modification 
factors for CNM [|B], as shown in Figure where they are compared to similar 
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Fig. 2. Nuclear dependence of J/ip produc- 
tion for three different energies vs X2 and Xf- 
Where xf = x\—X2 and X\ and xi are the mo- 
mentum fractions in d and Am respectively, a 
is a representation of the nuclear dependence 
in terms of a power law, i.e. a a = <J^A a . 
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Fig. 3. Nuclear modification factor 
versus rapidity for d + Au collisions. 
The yellow band bordered by black 
lines represents a fit to a model that 
contains EKS [ 7 shadowing and a dis- 
sociation cross section. 



data at lower energies. The lack of scaling with X2 shown in the left panel of the 
figure suggests that traditional shadowing models, which should have a universal 
X2 dependence, are not the dominant physics. The approximate scaling with xf 
(right panel) , at least for the lower energy data that extends to large xp , hints that 
initial-state energy loss or gluon saturation may be the dominant physics. 

In Figure [3] an approximate constraint using a simple CNM model (with shad- 
owing and dissociation) [ [10] is shown. This model can then be used to give an 
extrapolated constraint for Au + Au collisions, as shown in Figures 0] and [5] Clearly 
the d + Au data from 2003 used to constrain the CNM extrapolation here suffers 
from large uncertainties, and results in a large uncertainty for Au + Au collisions. 
For Au + Au at mid-rapidity the CNM band is almost consistent with the observed 
suppression - except for the most central collisions (n part ~ 340); while at forward 
rapidity the suppression seen for Au + Au is substantially stronger. The just com- 
pleted 2008 d + Au run has approximately 30 times more J/tp's than before and, 
once analyzed, will dramatically improve the knowledge of the CNM baseline in 
A + A collisions, and allow precision studies of the additional physics beyond CNM 
that comes from the hot-dense matter created in heavy- ion collisions. The CNM 
constraint is expected to narrow by approximately a factor of three with the new 
d + Au data. 
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Fig. 4. Extrapolation of the simple CNM Fig. 5. Extrapolation of the simple CNM 

model shown in Figure [3] to Au+Au colli- model shown in Figure [3] to Au+Au col- 
sions at mid rapidity. Results for both lisions at forward rapidity. Results for 
EKS[ [7j and NDSG[ [8] shadowing are both EKS[ [7] and NDSG[ [8] shadowing 
shown. are shown. 

4. How does the QGP affect Quarkonia 

Quarkonia are thought to be a definitive probe of the QGP through the screening 
process in the deconfined colored mcdium[ 11 . Different quarkonia states, because 
of their different binding energies, are expected to "melt" at different temperatures 
of the medium. E.g. in some lattice calculations the J ftp would melt at l.2Tc, 
but the T only at over 2Tc- Nuclear modification factors observed by PHENIX in 
Au + Au collisions are shown in Figure [6] The suppression at mid-rapidity is about 
the same as that observed for lower energies at the SPS[ 12J, despite the expectation 
that the hotter medium created at RHIC would cause a larger suppression. The 
suppression at forward rapidity is stronger than that at mid rapidity, and the ratio 
of the nuclear modification factors, forward/mid, shown in the bottom panel of the 
figure, reaches an approximately constant level of 0.6 for n part > 100. 

Several scenarios can be considered in trying to understand the observed trends: 
1) CNM effects, as discussed above, should always be accounted for as a baseline. 2) 
Sequential screening [[13] - where, as suggested by some lattice calculations, only the 
tp' and xc are screened and the J/ip itself is not - not at RHIC or at SPS energies. 
Then the observed suppression beyond CNM comes only from loss of the feeddown 
(~ 40%) from the two higher mass quarkonia states. 3) Regeneration models [[14]. 
where the large density of charm quarks created in the collisions (~ 20 in a central 
Au + Au collision) can produce charmonia in the latter stages of the expansion. 

In the sequential screening picture, if the CNM suppression at mid rapidity and 
the "melting" of the higher mass charmonia states was the same at RHIC and at 
the SPS, this would provide a natural explanation for the nearly identical suppres- 
sion at RHIC and the SPS. It would also agree with some lattice calculations that 
indicate no melting of the J/ip until over 2Tc [ US] . The stronger forward rapidity 
suppression seen at RHIC could then be explained by gluon saturation that gives 
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stronger forward suppression than that from standard shadowing models. For tra- 
ditional shadowing models the shadowing of the gluon from one nucleus is largely 
canceled by the anti-shadowing from the gluon from the other nucleus - resulting in 
an approximately flat rapidity dependence. For gluon saturation a " shadowing-like" 
effect is produced for the gluon in the small-x region, but no anti-shadowing for the 
other gluon, resulting in a stronger suppression at forward rapidity. Since screen- 
ing and gluon saturation might have different centrality dependences, it is unclear 
whether they would balance to produce the approximately flat ratio observed for 
n part > 100 (Figure^. 

An alternative is the regeneration picture, where the dissocation by the QGP 
at mid and forward rapidity would be similar, but the weaker suppression at mid 
rapidity would be due to regeneration effects being stonger here where the charm 
density is largest. In this case it would be an " accidental" compensation of screening 
and regeneration that leads to the same mid-rapidity suppression at RHIC and the 
SPS. At forward rapidity, where the charm density is smaller, the regeneration is 
reduced and stronger screening results. Again, whether the saturation in the forward 
to mid rapidity suppression could be reproduced by these two compensating effects 
is unclear. 

The regeneration mechanism depends on the square of the open-charm cross 
section, so it is critical to resolve the present uncertainties there. [ [To] Also, since 
charm has been shown to exhibit flow for moderate px values, one would expect 
J/tps that are produced by regeneration to inherit this flow. A first measurement 
of the J /ip> flow at mid rapidity is shown from part of the 2007 Au + Au data in 
Figure [7J but is clearly quite challenging, and so far is consistent with zero flow. 
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Fig. 6. Nuclear modification factor for 
Au+Au collisions at mid rapidity (red 
circles), and at forward rapidity (blue 
squares) versus centrality (top panel). In 
the bottom panel the ratio of the foward 
over mid rapidity nuclear modification 
factors from the upper panel is shown. 



"° ,3 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

Fig. 7. Flow of J ftps at mid rapFdffy V vs 
Pt (preliminary result from 42% of the 
2007 data), compared to several theoreti- 
cal models. 
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5. Summary and Future 

The suppression of J/ip production in Au + Au collisions at RHIC for mid rapidity 
is similar to that at lower energies, while for foward rapidity the RHIC suppression 
is stronger. Better cold nuclear matter constraints from the new d + Au data are 
needed to establish an accurate baseline and allow quantitative analysis of the QGP 
effects. Two theoretical pictures, 1) sequential suppression with gluon saturation 
and 2) dissociation and regeneration, appear to offer explanations of the observed 
trends. Higher luminosities and silicon vertex upgrades will enable much more 
quantitative studies in the next few years. Over 100,000 J/ips and 600 Ts are 
expected in a year with higher luminosities enabled by accelerator advances, while 
new silicon vertex detectors will allow explicit indentification of open-heavy and 
will improve both the background and mass resolution for the quankonia states - 
especially important to separate the ip' from the J/ tp at forward rapidity. 
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